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Abstract. For a rotating dust with a 3-dimensional symmetry group all possible metric 
forms can be classified and, within each class, explicitly written out. This is made possible 
by the formalism of Plebahski based on the Darboux theorem. In the resulting coordinates, 
the Killing vector fields (if any exist) assume a special form. Each Killing vector field may 
be either spanned on the fields of velocity and rotation or linearly independent of them. By 
considering all such cases one arrives at the classification. With respect to the structures 
of the groups, this is just the Bianchi classification, but with all possible orientations of the 
orbits taken into account. In this paper, which is part 1 of a 3-part series, all solutions are 
considered for which two Killing fields are spanned on velocity and rotation. The solutions 
of Lanczos and Godel are identified as special cases, and their new invariant definitions are 
provided. In addition, a new invariant definition is given of the Ozsvath class III solution. 

1. Introduction and summary. 

The theorem of Darboux presented in sec. 2 allows one to introduce invariantly defined 
coordinates in which the velocity field of a fluid (not assumed to be perfect) acquires 
a "canonical" form. In this paper it is further assumed that the fluid moves with zero 
acceleration and nonzero rotation. These assumptions result in a simplification of the 
metric tensor and in limitations imposed on the Killing vectors, if any exist. Within this 
special class of coordinates, any single Killing field may also be reduced to a "canonical" 
form, a different one in the case when it is spanned on the vector fields of velocity u a and 
rotation w a , and a different one when it is linearly independent of u a and w a . This gives 
rise to a classification of possible symmetries in rotating matter. 

When there exist three linearly independent Killing fields, the classification described 
above gives rise to a complete classification of all possible metric forms. With respect to 
the algebras of the symmetry groups, this is just the Bianchi classification, but with all 
orientations of the orbits in the spacetime taken into account. As a by-product, all the 
Bianchi types that are compatible with rotating dust, with any orientation of the orbit, 
are identified. 

In every case that emerges, the commutation relations of the algebra have been solved, 
resulting in explicit formulae for the Killing fields, and then the Killing equations have been 
solved, resulting in the formulae for the metric tensors compatible with the symmetry 
groups considered. The degree of success in solving the Einstein equations varied very 
strongly from case to case. In most cases, no headway was made. In some cases the Einstein 
equations have been integrated either to an autonomous set of first order equations or to 
a single nonlinear differential equation of second or third order. In a few cases solutions 
known earlier were identified in the present scheme and new invariant definitions for some 
of them were provided (those by Lanczos [1] and Godel [2] will be mentioned in this paper). 
In just one case a new solution was found. 



Since the number of cases is rather large, the results will be presented in three papers. 
The present paper 1 deals with the simplest situation when two of the Killing fields are 
spanned on velocity and rotation (the case of all three Killing fields being spanned on u 
and w is trivial - see sec. 6). 

The Darboux theorem was first applied as a tool for investigating the equations of 
motion and the Einstein equations by Plebahski [3]. He showed that if a perfect fluid is 
rotating and isentropic while the particle number is conserved, then a similar consideration 
to the one presented here applies. The approach of Plebahski was used by this author [4 
- 8] to find a large collection of stationary, cylindrically symmetric solutions of Einstein's 
equations. 

For a perfect fluid the assumptions of geodesic motion and nonzero rotation imply that 
the pressure is constant (see ref. 9). Hence, from the point of view of thermodynamics, 
geodesic and rotating perfect fluids are isentropic and fall within the class considered 
by Plebahski [3]. However, the approach based on the Darboux theorem applies to any 
timelike congruence that is of class C 1 and has zero acceleration and nonzero rotation. In 
particular, the velocity field of a rotating charged dust with zero Lorentz force, that was 
considered in several papers, has this property. The papers that discussed such solutions 
will be mentioned at the end of sec. 8; they are all within the same class of the classification 
introduced here. 

In sec. 2 the Darboux theorem is introduced. In sec. 3, the classification of first- 
order differential forms based on the Darboux theorem is applied to geodesic vector fields 
with rotation. When the vector field is the velocity field of a fluid, a class of preferred 
coordinates results, which shall be termed "Plebahski coordinates". In sec. 4, by way of 
an example, the consideration of sec. 3 is applied to the rotating dust solution of Stephani 
[10]. In sec. 5 it is shown that each Killing vector field that might possibly exist in a 
rotating dust spacetime is determined by two functions of two variables. If the Killing field 
is not spanned on velocity and rotation, then the Plebahski coordinates may be adapted 
to it so that it acquires the unique form k a = Sf . The Godel solution is used to illustrate 
the various forms of the Killing fields that may arise. 

In section 6, the consideration of sees. 3 and 5 is applied to the situation when there 
exist three Killing vector fields. When all three of them are spanned on u a and w a , the 
result is trivial: the group becomes two-dimensional, and this case is not considered here. 
When two of them are spanned on u a and w a while the third one is not, two cases arise 
that correspond to different Bianchi types (II and I) of the groups. These are investigated 
in sections 7 and 8. The solutions of Lanczos [1] and of Godel [2] emerge as special cases 
in both types. The Bianchi type II metrics are defined by a single third-order differential 
equation, the Bianchi type I metrics are determined by a set of autonomous first-order 
differential equations. 

Finally, in sec. 9, other invariant definitions are given: for the class III solution by 
Ozsvath [11] and for the solution of Godel [2]. The former results from the assumptions: 

1. The source in the Einstein equations is a geodesic, rotating perfect fluid. 

2. The rotation vector field is covariantly constant. 

The Godel solution, which is a subcase of this, emerges when it is assumed in addition 
that the shear of the fluid flow is zero. 

So far, Bianchi-type solutions of Einstein's equations with rotating source have been 
searched for and found by trial and error (often with nonperfect fluid sources, e.g. with 
viscosity or heat-flow). The results of the present series of papers show that, in the case 
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of dust source at least, the number of allowed possibilities is limited. It is hoped that the 
results will direct further research toward better-defined targets. 

2. The classification of differential forms of first order and the Darboux 
theorem. 

The Darboux theorem presented below exploits the simple fact that if a differential 
form q of first order is defined on an n- dimensional manifold M n , then its domain is not 
necessarily n- dimensional. Two cases are well-known: 

1. If q is a perfect differential so that q = df, where / is a scalar function on M n , then 
/ can be chosen as one of the coordinates, and the form becomes one-dimensional. 

2. If q has an integrating factor so that q = gdf, where / and g are independent scalar 
functions on M n , then / and g can be chosen as two of the coordinates, and then the 
domain of q is the two-dimensional (/, g)-surface. 

The Darboux theorem summarizes all the cases that can occur. It is based on the 
following classification (see also Ref. 3): 
Definition 

Let q be a differential form of first order. 

If Q21 '■= dq A . . . A dq [multiplied I times] / 0, but q A Q21 = 0, then q is said to be of 
class 21. 

If Q21+1 '■= q A Q21 / 0, but dQ2i+i = dq A Q21 = 0, then q is said to be of class (21 + 1). 

□ 

Then the following holds: 
The theorem of Darboux 

The form q is of class 21 if and only if there exists a set of 21 independent functions 

(€i,---,€i,vu---,vi) such that: 

q = md£i + m<%2 + ■■■■ + mdZi- (2.1) 

The form q is of class (21 + 1) if and only if there exists a set of (21 + 1) independent 
functions (r, £i, . . . , £j, rji, . . . , r)i) such that: 

q = dr + T&dfr + 7&df 2 + •••• + m^i- ( 2 - 2 ) 

A proof of this theorem can be found in Ref. 12. 

Evidently, the class of q cannot be larger than the dimension of the manifold on which 
q is defined. Hence, the Darboux theorem implies that in a four-dimensional spacetime V4 
any differential form of first order can be represented as: 

q = adr + rjd£, (2.3) 

where a, r, r\ and £ are scalar functions on V4 . 

Any vector field u a on V4 defines the following differential form: 

q u := u a dx a . (2.4) 

According to (2.3), in the most general case there exist scalar functions a,T,r] and £ such 
that: 

u a = ar, a +r]£,a ■ (2.5) 
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Note that the functions in (2.5) are not uniquely defined. Since we shall not use (2.5) in 
the most general case, we shall determine the nonuniqueness only in the subcase that is of 
direct interest to us (see sec. 3). 

For the most general case of (2.5), the four functions are independent in the sense that 
the Jacobian: 

9 ^ T ^0 (2.6) 



d(x° , x 1 , x 2 , x 3 ) 

Hence, they can be chosen as coordinates in the spacetime. In Refs. 4 and 7 it was shown 
that if u a is the velocity field of an isentropic perfect fluid in which the particle number is 
conserved, then a = 1/H, where H is the enthalpy per one particle of the fluid, and further 
limitations on u a follow from the particle number conservation. No other applications of 
(2.5) in the general case are known to this author. 

3. Geodesically moving fluids. 

To any timelike vector field u a normalized to unity (so that u a u a = 1) the formula 
from Refs. 13 and 14 may be applied: 

U a -/3 = U a U/3 + (Ta/3 + + (3.1) 

which gives rise to the well-known definitions of acceleration ii a , expansion 6, shear a a p 
and rotation u a p. In the signature (H ) used here, the projection tensor h a p is: 

h af 3 = g a p - u a up. (3.2) 
The following properties of ii a , a a p and u a p will be useful in further considerations: 

u a u a = 0, cr aj 3U^ = uJapuP = 0. (3.3) 

We shall assume from now on that u a is the velocity field of a fluid and that u a = 0, i. e. 
that the particles of the fluid move on geodesies. Then, from (2.5) we have: 

w «/3 = ^( u a,/3 ~ U/3,a) = T ^ ~<7,a T,p +rj,/3 £,a ~V,a £>/3 ), (3.4) 

and from (3.3) we have: 

{vPa,p )t, q -{u P T,f3 )a, a +(n /3 ??, /3 )^, Q )rj, a = 0. (3.5) 

There are two possibilities now: 

I. At least one of the four scalar products in (3.5) is nonzero. In this case (3.5) implies 
that at most three of the functions (a, t, 77, £) are independent, and so the form (2.4) will 
not be of class 4. 

II. All the four scalar products are zero. However, this means that the gradients of 
(<r, t, £,77) are all confined to the 3-space orthogonal to u a , i.e. that there is a functional 
relation among these four functions. Again, the form (2.4) cannot be of class 4. 

Hence, for a geodesically moving fluid the form (2.4) is of class at most 3, i.e. at most 
3 independent functions r, rj, £ exist such that: 

U a = T, a +??£, Q . (3.6) 
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From here on, the reasoning used in Refs. 3 and 4 applies almost unchanged. With (3.6) 
we have in (3.4): 

= ^(V,/3 £,a ~V,a £,p), (3.7) 

and in (3.5): 

(u\ ( 3)^ a -(u^,p)r ] , a =0. (3.8) 

There are again two possibilities: 

I. Either (u@r],Q ) and (u^^,p ) do not vanish simultaneously, and then (3.8) implies that 
r\ and £ are functionally related, in which case (3.6) implies that u a is a gradient of a 
function, and so u a p = 0. 

II. Or £ and rj are not functionally related, in which case: 

u^, p = vPr)^ = 0, (3.9) 

and Ld a /3 0. We shall be interested only in the second case. The functions {r, £, 77} in 
(3.6) are determined up to the following transformations: 

t = Fg,rf), r, = G{?,rf), r = r' - S{g ,r/), (3.10) 
where the functions F and G must obey the equation: 

F,£> G, v i —F, v i G,£ = 1, (3-11) 

and then S is determined by: 

= GF,£> —rj', S, v i = GF, V > . (3-12) 

Eq. (3.11) is the integrability condition of eqs. (3.12) and it ensures that the Jacobian 
of the transformation (3.10) equals 1. It follows that one of the functions {F, G} can be 
chosen arbitrarily, the other one is then determined by (3.11) and S is fixed up to an 
additive constant by (3.12). The inverse transformation to (3.10) is of exactly the same 
form, with the corresponding F, G and S obeying (3.11) and (3.12). 

Let us now make the additional assumption that the number of particles of the fluid is 
conserved, i.e.: 

(v^m/% = 0, (3.13) 

where g is the determinant of the metric tensor and n is the particle number density. This 
equation is a necessary and sufficient condition for the existence of a function Q such that: 

^nu a = e a ^ s ^ Vn C, 5 . (3.14) 
Note that (3.6) and (3.9) imply that: 

u a r, a =l, (3.15) 

and then eq. (3.14) implies that: 
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e a ^ S r, a £ )/3 C,S = g^t^) =^n^0. (3.16) 
Eq. (3.14) implies also that: 

u a (, a = 0. (3.17) 
The function £ is determined by (3.14) up to the transformations: 

C = e + T(£,rf), (3.18) 

where T is an arbitrary function. Eq. (3.16) certifies that {r, £, r/, can be used as coordi- 
nates in the spacetime. If they are chosen as the {x . X . X • iE } = {t, x, y, z} coordinates, 
respectively, then eq. (3.6) implies: 

u = 1, u\ = y, u 2 = u 3 = 0. (3.19) 

We will use these coordinates throughout the remaining part of the paper and call them 
"Plebahski coordinates". Eq. (3.16) implies now: 

g = -n~ 2 , (3.20) 

and eq. (3.14) implies: 

u a = SS, (3.21) 
i.e. the Plebahski coordinates are comoving. The rotation vector defined by: 

w a = -(l/^) £ a ^ s u^s (3.22) 

assumes the form: 

w a = nd$. (3.23) 

Eqs. (3.19) and (3.21) imply that: 

5oo = 1, 5oi = V, 902 = 903 = 0, (3.24) 
and also that the only nonvanishing components of the rotation tensor are: 

uj 12 = -u 2 i = 1/2. (3.25) 

Note that, in contrast to Refs. 3 and 4 where barotropic perfect fluids were considered, we 
have not assumed anything about the form of the energy-momentum tensor so far. 

If we now assume that the fluid is perfect, then we conclude from the equations of 
motion T a ^-p = that either oj = or p = const (see also ref. 9). This means that 
a geodesic perfect fluid can be rotating only if it is in fact dust; the constant p can be 
reinterpreted as the cosmological constant. In this case, the energy-density obeys the 
conservation equation (y/— geu a ) :Ct = and eq. (3.13) need not be assumed separately^. 

For dust, results closely analogous to (3.19) - (3.25) were obtained by Ellis [15], by adapting an or- 
thonormal vector basis and a coordinate system to u a and w a . Of the exact solutions with nonzero rotation 
found by Ellis most, but not all, do belong to the collection considered in this series of papers. They will 
be described in paper 2 
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4. Example: the Stephani solution [10]. 

The Stephani metric with p = const 7^ (eq. (4.22) in Ref. 10) is not in fact a perfect 
fluid solution, as was found out while trying to construct the Plebahski coordinates for itQ. 
Therefore, we shall consider only the dust solution, eq. (4.8) in Ref. 10. In the original 
notation except for the signature, the solution is: 

ds 2 = r, ab dx a dx h - ^{dx 1 ) 2 , (4.1) 
where a, b = 0, 2, 3, r) ab = diag (1, -1, -1) and: 

N:= ^MlnT + g a x a + h, 

T 2 :=rj ab (x a -r)(x b -f b ), (4.2) 

the functions M(x 1 ), f a (x 1 ), g a (x 1 ) and /i(x x ) all being arbitrary. The velocity field and 
the energy-density of the dust are, respectively: 

u a = T, a , m = 0, (4.3) 
Ke = M/(NT 2 ). (4.4) 
The formula for the velocity field can be writen as: 

u a = T, a -T, 1 x 1 , a , (4.5) 
which immediately suggests the choice of the Plebahski coordinates of sec. 3: 

t = T, Z = x l , r, = -T, 1 . (4.6) 
Eq. (3.14) defining £ is here: 

MT- 2 V ab T, b = -e alcd T, lc (, d . (4.7) 
The following identity is useful in calculations: 

T, 2 - T, 2 2 - T, 3 2 = 1. (4.8) 

Using this, one can verify that only two of the three equations (4.7) are independent. The 
solution of (4.7) is: 

( = M[U + \(T, X )], (4.9) 
where A(T,i ) is an arbitrary function, and U is defined by: 

U = J [f a (1 - X 2 ) - f° a Y{X) + f 2 a XY{X)}- l dX, (4.10) 
the function Y(X) being determined by: 

T,i (1 - X 2 - Y 2 fl 2 = -/% +/ 2 ,x X + f a Y. (4.11) 

2 The error is deeply hidden and so far could not be corrected. I am grateful to H. Stephani for cooperation 
on this point. 
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In the integral (4.10), the coordinate x 1 and the quantity T,\ are to be treated as parameters 
independent of X. After the integral is calculated, one should substitute for T,i from (4.11), 
while X and Y are to be replaced by: 

X = (x 2 - f 2 )/(x° - f), Y = (x 3 - f)/(x° - f), (4.12) 

The integral in (4.10) is expressible in terms of elementary functions, but the result is 
rather complicated. 

As can be seen, the solution (4.1) - (4.4) becomes very complicated in the Plebahski 
coordinates and it is unlikely that it could be found in such a form from the Einstein 
equations. Thus, although the coordinates are invariantly defined, they are not necessarily 
convenient for explicit calculations. 

A collection of other solutions represented in the Plebahski coordinates can be found 
in the extended version of Ref. 6. 

5. The Killing vector fields compatible with rotation. 

We shall assume that the symmetries of the spacetime (if any exist) are inherited by 
the source, i.e. that if the Lie derivative of the metric tensor g a p along the vector field k a is 
zero, £k9a/3 = 0, then the velocity field and the particle number density are also invariant: 
£kU a = = £kn. (For a pure perfect fluid source the inheritance is guaranteed.) It follows 
that the rotation tensor must also be invariant, £k^a/3 = 0. 

In consequence of (3.21) the equation £kU a = k^u a ,^ —u^k a ^ = implies that: 

k a , t = 0. (5.1) 

In consequence of (3.23) and of the assumption £kn = k a n, a = 0, the equation £kW a = 
implies: 

k a , z =0. (5.2) 
The equation £k^ a (3 = 0, in consequence of (3.25) implies: 

k ix +^ iy = 0) (5-3) 

and the equation £ku a = 0, in consequence of (3.19) implies: 

k° , x = —k 2 — yk} , x , k°,y = —yk},y . (5-4) 

(The equations £kU a = and £kU a = provide independent pieces of information because 
the equations £k9a/3 = have not been used so far). Eq. (5.3) is the integrability condition 
of eqs. (5.4). The general solution of eqs. (5.1) - (5.4) is: 

k° = C + <j) - y<p, y , k l = (j), y , k 2 = -(j), x , k 3 = X, (5.5) 

where (p(x,y) and X(x, y) are arbitrary functions and C is an arbitrary constant. Sym- 
metries need not be present, in fact the Stephani [10] solution considered in sec. 4 is an 
example of a rotating dust solution with no symmetries. In this case <p = X = C = 0. 
However, if any symmetries are present, then the Killing vector fields must have the form 
(5-5). 

Suppose that <fi is not a constant, i.e. that a Killing vector field k a exists that has a 
nonzero component in the x- or y-direction (in invariant terms this means that the vector 
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field k a is not spanned on the vector fields of velocity, u a , and rotation, w a ). We can then, 
within the Plebahski class defined in sec. 3, adapt the coordinates to k a in such a way 
that k a = 5f , i.e. so that the metric becomes independent of x'. From (3.10) - (3.12) and 
(3.18) the transformation functions are: 

t' = t-S(x,y), x' = F(x,y), y' = G(x,y), z' = z + T(x,y), (5.6) 
where T is arbitrary, while F, G and S obey: 

F, x G, y —F, y G, x = 1, S, x = GF, X —y, S, y = GF, y . (5-7) 
In order to lead to k a ' = 5f' the functions F, G and T must obey in addition: 

-(<P + C) + GF, X fry -GF, y <f), x = 0, 

F, x 4>,y—F,y4i, x =l, G, x (j), y — G,y (j), X = 0, (5.8) 

T, x <fi,y -T, y </), x = -A. (5.9) 

The unique solution of eqs. (5.8) is G = 4> + C, which obeys (5.7) as well (in virtue of 
the second of (5.8)). Eq. (5.9) simply defines the accompanying T which is seen to exist 
always. Since <j) was assumed nonconstant, the transformation is nonsingular (in fact its 
Jacobian equals just 1), and results in 4> = y in the new coordinates. As already noticed, 
the metric becomes independent of x after the transformation. This property is preserved 
by the transformations (5.6), but with F, G, S and T restricted now by: 

G = y, F = x + H(y), T = T(y), S = J yH, y dy + A, (5.10) 

where A is an arbitrary constant and H, T are arbitrary functions. The functions given by 
(5.10) fulfil (5.7) identically. Note that the transformation to k a = 5f exists irrespectively 
of any possible functional relation among A, (p and C in (5.5); the only case when it fails is 
4> = const. 

A solution of the Einstein equations may have more than one Killing vector field of the 
form (5.5). In that case, the transformation (5.8) - (5.9) changes only one of them to the 
preferred form, the others will preserve their more complicated appearance. An example of 
this situation is the Godel solution [2] transformed to the Plebahski coordinates, see Refs. 
4 and 5: 

ds 2 = {dt + ydxf - ^y 2 dx 2 - (Key^dy 2 - 2Ke~ l dz 2 , (5.11) 
where k = 87rG/c 4 and e is the energy-density related to the cosmological constant by^: 

A=^Ke. (5.12) 

The symmetry group of this solution is 5-dimensional, the independent 1-parametric sub- 
groups were given in Ref. 4. Those connected with nonconstant <j) in (5.5) are the following 
three: 

3 Note that if the cosmological constant is reinterpreted as pressure, then the resulting perfect fluid has 
the equation of state e = p. Hence, the Godel solution may have been the first example considered in the 
literature of a "stiff perfect fluid", now familiar from the studies of solution-generating techniques, see e.g. 
Verdaguer [16]. 
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x = x' + Si, 
x = e~ S2 x', y = e S2 y\ 



(5.13) 
(5.14) 



t = t'+ (2V2/K) arctan[V2s 3 (Ky'y l (l - S32T 1 ], 

x = [K 2 x'y> 2 (l - s 3 x') - 2s 3 ]/[2a! + K 2 y' 2 (l - s 3 x') 2 ], 

y= (l-s 3 x') 2 y' + 2s 2 /(K 2 y'), (5.15) 

where s\,S2 and S3 are the group parameters and K := (Ke) 1//2 . The Killing vectors are, 
respectively: = 5f (corresponding to C = A = 0, 4> = y, the one constructed in (5.8) 
- (5.9)), k? 2 \ = —x5f + y&2 (corresponding to <j) = —xy) and k^ = A(K 2 y)~ 1 5 l Q + [x 2 — 
2/(Ky) 2 ]5f - 2xy5% (corresponding to <fi = x 2 y + 2/(K 2 y)). 

6. The algebra of three Killing fields. 

Suppose that three Killing vector fields exist and all three are spanned on u a and w a , 
so that (f) = const in (5.5) for each of them, i.e.: 

k( i} = C i 6S + \ i (x,y)6%, i = 1,2,3. (6.1) 

From the Killing equations one can then easily conclude that constants a\, «2 and a 3 exist 
such that aikn\ +02^(2) + a 3^(3) = 0, i.e. the symmetry group is in fact two-dimensional. 
Hence, no three- dimensional symmetry group with the generators (6.1) exists^; for a three- 
dimensional group at least one of the generators must be linearly independent of u a and 
w a at every point of the spacetime region under consideration. 

In sections 7 and 8 we shall consider the situation when exactly one generator, k?^, is 
everywhere linearly independent of u a and w a , while the other two, k%\ and k? s \, are of 
the form (6.1). In agreement with the result of sec. 5, the Plebahski coordinates can be 
adapted to k%\ so that: 

fcfl) = ft, (6-2) 

while: 

*$) = C ^o + A 2 (x,y)5 3 a , kf 3) = C 3 Sg + A 3 (x,i/)<5£, (6.3) 

and the coordinate transformations preserving (6.2) and (6.3) are (5.10). Note that C2 
and C3 cannot vanish simultaneously because otherwise the Killing equations immediately 
imply that either k? s \ = const k%\ (in which case the symmetry group is two-dimensional) 
or the metric is singular. However, with no loss of generality we can assume that: 

C 2 + = C 3 (6.4) 

because the Killing vector fields are determined up to linear combinations among them. 
Hence, if initially C2 / / C3, then we take k'fy = — (C 3 /C2)k^ instead of as 
the basis generator. If initially C2 = 7^ C3, then we exchange the labels "2" and "3". 

4 The algebra (6.1) corresponds to a three-dimensional group that has two-dimensional orbits, and it 
turns out that in the case considered the group has to be two-dimensional as well. As will follow from 
the whole of the present work, three-dimensional symmetry groups with two-dimensional orbits just do not 
exist for rotating dust. 
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We shall keep the choice (6.4) throughout sections 7 and 8. 

7. The generators, the Killing equations and the Einstein equations for a 
Bianchi type II algebra. 

All the cases that arise follow as limits at different stages of calculation from the generic 
case A2 / 0, and we shall consider it first. The commutators of the Killing vectors are: 

[k(ih k (2)] a = (AWA 3 )£f 3) , [k ( 2),k {3) ] a = 0, 

[fc(i),*(3)] a = (A 3) x/A 3 )Af 3) . (7.1) 
The Killing vector fields will thus form a Lie algebra when: 

X 2 ,x = b\ 3 , X 3jX = cA 3 , (7.2) 

where 6 and c are arbitrary constants. The case c / has to be considered separately. 
Then: 

A3 = P(y)e cx , A 2 = {b/c)f3{y)e cx + a(y), (7.3) 

where a(y) and (3{y) are arbitrary functions. However, in this case we can take = 
^(2) ~~ (^/ c )^(3) as the new basis generator instead of fc^, and the result is equivalent to 
assuming 6 = 0. Hence, with c / 0, we can take 6 = with no loss of generality. 
The further procedure consists of the following steps: 

1. Adapt the coordinates to the Killing fields to make them as simple as possible. 

2. Solve the Killing equations for the components of the metric tensor. 

3. Simplify the metric as far as possible by coordinate transformations. 

4. Solve the Einstein equations. 

The coordinate transformations in steps 1 and 3 in general lead out of the Plebahski 
class defined in sec. 3. 

This procedure will be presented in some detail below. In the present case the result 
from the Einstein equations is: either c = or there is no rotation. Since we are interested 
in rotating solutions only, the case c / = 6 need not be followed further. We thus assume: 

c = 0. (7.4) 

Then: 

As = /%), X 2 = bP(y)x + a(y). (7.5) 

The algebra of the Killing vector fields is of Bianchi type II when 6/0 and of Bianchi 
type I when 6 = 0. 

In order to simplify the Killing vectors we now transform the coordinates as follows: 

(t',x',y') = (t,x,y), z> = -( a /C 2 )t + z/p. (7.6) 

The transformation is not of the form (5.10), so the new coordinates do not belong to the 
Plebahski class, and the forms of velocity, rotation and the metric will no longer agree with 
(3.19) - (3.25). The Killing vector fields in the new coordinates become: 
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fcfi) = k (2) = *o + bx5%, £f 3) = 5%, (7.7) 

while the velocity and rotation fields become: 

u a = 6g - (a/C 2 )5%, w a = (n//3)8§. (7.8) 

The transformed metric is independent of x and z, while the Killing equations for 
imply: 

5oo = 1 + (a/C 2 ) 2 h 33 (y), g 01 = y + (a/C 2 )g 13 , 

902 = 0, c/03 = {a/C 2 )h 33 (y), 

9ii = h 33 {y){bt) 2 - 2h 13 (y)bt + h n (y), 

912 = ~h 23 (y)bt + h 12 (y), g 13 = -h 33 (y)bt + h 13 (y), 

922 = h 22 (y), g 23 = h 23 (y), g 33 = h 33 (y), (7.9) 

where hij(y),i,j = 1,2,3, are arbitrary functions of y, to be found from the Einstein 
equations. 

The orbits of the symmetry group are now the hypersurfaces y = const. In order to 
follow the standard technique of the Bianchi-type spaces we should now carry out a coordi- 
nate transformation that preserves (7.7) and makes the y-coordinate curves orthogonal to 
the group orbits, so that g'^ = g[ 2 = g 23 = after the transformation. This step is not in 
fact necessary for solving the Einstein equations (in general it only reshuffles the unknown 
functions without elliminating any of them), but in the case under consideration it leads 
to a simplification. The transformation is: 

t = t' + f (y'), x = x' + f 1 (y'), y = f 2 (y'), 

z = z' + bh{y')t' + h{y'), (7.10) 

where f a (y') must obey: 

fo,y' = —hh,y' + ( a /C 2 )h 23 f 2j y/, 
(-/I + hu - hl 3 /h 33 )f hy > + (h 12 - h 13 h 23 /h 33 + af 2 h 23 /C 2 )f 2 , y > = 0, 

h, y ' = -{a/C 2 )f ,y> + (bfo - h 13 /h 33 )f hy > - {h 23 / h 33 ) f 2tV > . (7.11) 

The equations are well-defined because: 

I- ^33 / 0; otherwise the rotation vector would be null, which is a physical impossibility. 

II. —/I + fin — ^13/^33 / 0; otherwise the determinant of the metric tensor becomes 
positive; i.e. the metric acquires an unphysical signature. 

Eqs. (7.11) are to be understood as follows. The function f 2 (y') can be chosen ar- 
bitrarily, therefore we choose it so that g' 22 = —1 after the transformation. With f 2 (y') 
thus chosen, fi(y') is found from the second of (7.11), then fo(y') is found from the first 
of (7.11), and finally f 3 (y') is found from the third of (7.11). 

After the transformation the metric becomes (primes dropped): 

5 oo = l + (a/C 2 + bfi) 2 h 33 , 
901 = Y(y) + (a/C 2 + bfi){-h 33 bt + h 13 - bf h 33 ), 
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502 = 512 = 523 = 0, 

503 = (a/C 2 + 6/i)/i33, 
511 = h 33 b 2 (t + / ) 2 - 2h 13 b(t + /„) + /in, 
513 = ~h 33 b(t + /o) + /ii3, 

522 = -1, 533 = h 33 (y), (7.12) 

where hij(y),i,j = 1, 2, 3, / (y), /i(y), Y(y) = f 2 (y) and a(y) are functions to be found 
from the Einstein equations, and b and C 2 are arbitrary constants, C 2 ^ 0. 

For convenience in calculations we introduce the new functions G(y), A(y), k\ 3 (y) and 
F(y) by: 

533 = -G\y), a/C 2 + bf x = A(y), h 13 = -G 2 (k 13 + bf ), 

h n =Y 2 -F 2 - k 2 3 G 2 + b 2 G 2 fl - 2bf k 13 G 2 . (7.13) 
The velocity field in the coordinates of (7.12) - (7.13) is: 

u a = 6%-A5%. (7.14) 

Now the metric form is: 

ds 2 = (dt + Ydx) 2 - {Fdxf - dy 2 - G 2 [Adt - (bt - k 13 )dx + dz] 2 . (7.15) 

The components of the Einstein tensor will be referred to the orthonormal tetrad of forms 
e l = e l a dx a ,i = 0,1,2,3, uniquely implied by (7.15). Note that e° = u a dx a . Hence, the 
Einstein equations are: 

Goo = (87rG/c 4 )e, 

Gn = G 22 = G 33 = A, Gn = when i + j, (7.16) 

where e is the energy-density and A is the cosmological constant. 

The equation G\ 2 = implies that bA, y = 0. The case 6 = will be considered sepa- 
rately below, so we take here: 

A = const. (7.17) 

Then G02 = implies: 

k 13 = const. (7.18) 
We can then carry out the coordinate transformation: 

z = z' -At- k 13 x, (t, x, y) = (t', x', y'), (7.19) 
which has the same result as if: 

A = fci 3 = 0, (7.20) 

and we shall assume this from now on. The metric is still independent of x and of z, while 
A = k± 3 = implies g± 3 = 0, i.e. the Killing vectors fc^ = 5f and k^ = 5% are orthogonal 
to each other. The equation Gqi = then has the integral: 
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Y, y G/F = B= const, (7.21) 

and we can assume B / because rotation would be zero with B = = Y, y . 

At this point, only the diagonal components of the Einstein tensor are still nonzero, of 
which Goo just defines the energy- density, and the other three are functionally dependent 
(i.e. if Gn = A = G22 are fulfilled, then so is G33 = A). They determine F(y) and G(y). 

It is convenient to introduce Y(y) as the new variable. The equation Gn + G22 = 2A 
can then be written, with the help of (7.21), as: 

(F 2 G,y/G), y = 2AG 2 /B 2 -^, (7.22) 

and so: 

F 2 = (C-^ + 2 A| G 2 dY)G/G, Y , (7.23) 

where C is a new arbitrary constant (we can assume G,y / because G,y = implies 6 = 
from Gn — G22 = 0, and 6 = will be considered separately). Using (7.23) in G22 = A we 
obtain the following integro-differential equation that determines G: 



-^6 2 GG,y + I(5/G) 2 (G - \y + 2 A J G 2 dY) 2 (G, Y /G - G, YY /G, Y ) = 0. (7.24) 

In the special case A = this becomes an ordinary second-order differential equation. 
It is easy to get rid of the integral by transforming (7.24) appropriately and differentiating 
the result by Y(in this way a third-order differential equation for G(Y) is obtained) or 
by introducing the new variable u(Y) by dY/du = 1/G 2 (this results in a second-order 
equation for G{u)). However, no progress toward solving (7.24) results in either case. 

With the help of the equations Gn = A = G22 the formula for energy-density may be 
simplified to: 

(8vrG/c 4 )e = (B/G) 2 - (6G) 2 - 2A. (7.25) 

Note that the solutions considered here have a meaningful limit 6 = 0. 

When G = const, eqs. (7.23) and (7.24) no longer apply and one has to go back to the 
Einstein equations. They imply: 

G 2 = B 2 /{AA) (7.26) 

(i.e. necessarily A > 0) and: 

F 2 = ^Y 2 + DY + E, (7.27) 

where D and E are constants. If Y is chosen as the new coordinate in place of y, then 
from (7.21) and (7.26) the metric component gyy is: 

9YY = —(G/BF) 2 = -1/(4AF 2 ) (7.28) 

and the resulting metric is the Godel solution (see Ref. 4). Note that G = const is 
equivalent to e = const, see eq. (7.25). 
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When G,y ^ = b, eq. (7.24) implies G = e DY+E , and this leads to the Lanczos 
solution (see Ref. 4). 

These derivations of the Lanczos and Godel solutions lead to their invariant definitions 
that are based on weaker assumptions than the definitions known so far: 

1. The source in the Einstein equations is a rotating dust. 

2. The spacetime has a 3-dimensional symmetry group. 

3. Two of the symmetry generators are spanned on the vector fields of velocity u a and 
rotation w a , while the third one is linearly independent of u a and w a at every point. 

4. The generators form a Bianchi type II algebra. 

5. In the solutions of the Einstein equations, the Bianchi type I limit is taken of the 
Bianchi type II symmetry. 

6. The Godel solution results when the matter-density is constant, the Lanczos solution 
results when the density is not constant. 

The generalization with respect to the earlier definition is contained in point 3: in 
previous derivations the two generators were assumed to be collinear with u a and w a , 
respectively, from the beginning. 

8. The generators, the Killing equations and the Einstein equations for a 
Bianchi type I algebra. 

We shall consider the case b = c = in (7.1) - (7.2). The reasoning up to eq. (7.16) 
applies also here, but (7.17) no longer follows. Instead, the equation G13 = can be 
integrated with the result: 

ki3,y = BF/G 3 — YA, y , (8.1) 
where B is an arbitrary constant; the equation Got = can be integrated to: 

Y, y =(C -BA)F/G, (8.2) 
where C is an arbitrary constant; and the equation G03 = can be integrated to: 

A, y =(BY -D)/(FG 3 ), (8.3) 

where D is one more arbitrary constant. 

At this point, only the diagonal components of the Einstein tensor survive, and Goo = 
(87rG/c 4 )e — A just defines the energy- density. The equations Gn = A = G22 = G33 can 
be written as[] : 

B , C-BA„ „ 2BY-D . 

~~4FG ' y + ~~AFG~ y /G ~ TFG = A ' (84) 

B k C ~ BA Y 1 Em G >y D A -A rstt 

3D 

4pg**» ~ -JFG- Y >y +F >yy > F + ^fg Ay = A ' (8 ' 6) 

The set (8.4) - (8.6) can be integrated to a first order set. Subtracting (8.6) from (8.4) and 
multiplying the result by FG we obtain an equation that is easily integrated to: 



5 In order to arrive at this form, one has to calculate B from (8.1) and replace one factor B in B 2 by the 
resulting expression; then replace one Y, y in Y 2 from (8.2) and replace one A, y in A 2 from (8.3). 
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FG, y —GF, y -Bk 13 - - BAY + -CY - -DA = E= const. (8.7) 

Now adding (8.3) and (8.4), and multiplying the result by FG we obtain another integrable 
equation whose integral can be written in the form: 

FG, y = ^Bk 13 + X - BAY - X -CY + 2\j FGdy + H , (8.8) 

where Hq is an arbitrary constant. The integral can be calculated if the new variable u(y) 
is introduced by: 

dy/du = 1/{FG). (8.9) 

From (8.7) and (8.8) it follows that: 

GF, y = -^Bk 13 - \-DA -E + 2A f FGdy + Hq. (8.10) 

In the set (8.4) - (8.6) there remains one equation that has still not been used. However, 
at this point it merely introduces a relation between the arbitrary constants, i.e. implic- 
itly defines Hq in terms of the other constants. This is seen as follows: substitute for 
ki 3t y,Y, y , F, y , G ,y and A, y from (8.1) - (8.3), (8.8) and (8.10) in (8.5), thereby obtaining 
an algebraic equation (i.e. one without derivatives). Differentiate it by y and elliminate the 
derivatives in the same way again. What results is an identity = 0. Hence, the left-hand 
side of (8.5) is identically constant in virtue of the other equations. 

In terms of the variable u from (8.9), eqs. (8.1) - (8.3), (8.8) and (8.10) form an 
autonomous set of first-order equations that can be investigated further by qualitative 
methods (see e.g. Ref. 17). This is left as a subject for a separate study. 

In analogy with the Bianchi type I spatially homogeneous (nonrotating) dust solutions 
(see Ref. 18) one might expect further progress by adapting the coordinates suitably (in 
the case considered in Ref. 18, the metric can be diagonalized). However, this author was 
not able to achieve any such progress. 

The functions A(y) and k\ 3 {y) have invariant meaning: they are proportional to the 
scalar products of the Killing vectors (see eqs. (7.7) and (7.15) with b = 0): 

A = -g a pkf 2) k^/G 2 , k 13 = -g^kf^/G 2 (8.11) 

(note that G 2 = —g a/ 3k^k^, i.e. it is a scalar, too). Hence, A = and k\ 3 = are 
invariant properties. Note that ^4 = implies, through (8.3), that either Y = const (in 
which case there is no rotation) or B = D = 0. In the latter case, k± 3 = const and the 
coordinate transformation z = z' — k± 3 x leads to ki 3 = in the new coordinates. With 
A = k\ 3 = 0, the Lanczos and Godel models result from the Einstein equations as the 
only solutions. Hence, one more invariant definition of these models follows, similar to the 
six-point definition at the end of section 7. Points 1, 2, 3 and 6 remain unchanged, while 
points 4 and 5 are replaced by: 

4'. The generators form a Bianchi type I algebra. 

5'. From the two generators spanned on u a and w a , two linear combinations can be 
constructed that are orthogonal to each other. 

Point 5' is equivalent to the existence of coordinates in which A = 0. 
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Note that the Bianchi type I models considered in this section are more general than 
the Bianchi type I limit of the models from sec. 7; those from sec. 7 had A = ki% = in 
virtue of Einstein's equations. 

The assumption k%s = (i.e. Sapk^kfL = 0) alone does not lead to any immediate 
progress in solving the Einstein equations. 

The Lanczos solution was originally derived in Ref. 1 (an English translation, Ref. 19, 
is now available), and rediscovered in Ref. 20. Its limit of zero cosmological constant was 
rediscovered in Ref. 21 as the cylindrically symmetric subcase of a family of stationary 
axially symmetric solutions. Geometrical and physical properties of the Lanczos solution 
were discussed in Ref. 1, and, in a more modern language, also in Ref. 22 (the latter only 
for the case A = 0). 

Coordinate transforms of the Godel solution were published as new solutions in Refs. 
23 and 24 (concerning Ref. 23 see also Ref. 25). 

A metric form that is a modest generalization of the Godel solution (it has two unknown 
functions of one variable in place of Godel's e x and e 2x ) came to be known as "Godel-type 
metric" and became the subject of a rather large number of papers; the activity seems to 
have started with Ref. 26, one of the most recent appearances of it is Ref. 27. However, it 
was proven already in Ref. 28 that the only perfect fluid solution with this metric is the 
Godel solution itself; indeed, all other "Godel-type solutions" have various nonperfect fluid 
sources, and therefore they do not show up in the scheme considered here. 

As mentioned in sec. 1, several authors considered rotating charged dust solutions under 
the additional assumption that the electromagnetic field exerts no force on the charged 
dust particles, i.e. that F^u 11 = 0. These solutions were all derived with another, rather 
natural assumption: that all charges are attached to dust particles so that no currents are 
present apart from the one created by the dust flow. Those solutions are found in Refs. 
29 - 36. The one in Ref. 29 has only two-dimensional symmetry, so it could not come 
up in this investigation. The remaining ones are stationary and cylindrically symmetric 
and would have shown up here, had we allowed charges and electromagnetic fields in the 
source. They have the following properties: 

The one from Ref. 30 becomes a vacuum solution in the limit F^ = 0. 

The one from Ref. 31 does not allow this limit at all. 

The limit F^ = of the solution from Ref. 32 is the Minkowski metric. 

The Som - Raychaudhuri solution [33] reproduces the A = subcase of the Lanczos 
solution when F^ u = 0. 

The first of the six solutions by Banerjee and Banerji [34] reduces to the Godel solution 
when Fn V = 0. The other five behave as follows: 2 and 6 become vacuum solutions when 
Fav = 0, no 5 becomes the Minkowski spacetime, no 3 does not allow this limit at all, and 
no 4 has has a two-dimensional symmetry group. 

Both solutions by Mitskievic and Tsalakou [35] are generalizations of the Godel solution; 
the first one of them is in addition a generalization of the full (A ^ 0) Lanczos solution^]. 
In the limit A = 0, the first solution reduces to the one by Som and Raychaudhuri [33]. 

The two solutions from Ref. 36 are coordinate transforms of those from Ref. 35. 

Three other generalizations of the Godel solution exist in the literature that have zero 
acceleration. Two were provided by Raval and Vaidya [37]; the first of them is station- 
ary, the second expanding, both have anisotropic pressure. The third is the solution by 

6 In fact, the second solution has nonzero pressure gradient that remains nonzero even after the limit 
Fpv = is taken. Another limiting transition, given in the paper, reduces the solution to Godel's. 
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Reboucas [38] in which the source is a free electromagnetic field (see also Ref. 39). The 
metric of the Reboucas solution is the same as that in the first Banerjee - Banerji solu- 
tion. This coincidence was explained by Raychaudhuri and Guha Thakurta [40]: The two 
electromagnetic fields (one generated by a current, the other source- free) are related by a 
point-dependent duality rotation. 

9. Another invariant definition of the Godel and Ozsvath class III solutions. 

Assumptions about invariant properties of the velocity field of matter usually lead to 
progress in solving the Einstein equations; the most impressive example were the shearfree 
normal models of Barnes [41], where a large class of solutions resulted from the assumptions 
of zero shear and zero rotation in a perfect fluid source. Inspired by this, one can try to 
make assumptions about other vector fields characterizing fluid sources, e.g. the rotation. 
Indeed, it turns out that the assumption: 

w a;P = 0, (9.1) 

i.e. the rotation field being covariantly constant, together with the assumption of geodesic 
motion of a perfect fluid source, leads uniquely to two solutions of Einstein's equations. 
However, both of them were obtained before by other methods. One is the Ozsvath class 
III metric [11], originally identified as one of the solutions that are homogeneous in four 
dimensions; the other is the Godel solution [2] which is the shearfree limit of the Ozsvath 
solution. 

From (9.1) and from the Ricci identity 2w a .[^ = R p a p^w p one obtains for the Ricci 
tensor: 

R%w p = 0, (9.2) 
and then from the Einstein equations for a perfect fluid: 

G a/ 3 + Ag a/3 = n[(e + p)u a u p - pg a p], k = 8vrG/c 4 , (9.3) 
and from u a w a = one obtains: 

A = l«(e-p). (9.4) 

In the case A = 0, this is the well-known "stiff perfect fluid". Eq. (9.4) is a necessary 
condition for (9.1) when the source is a perfect fluid. 

As stated at the end of sec. 3, if the perfect fluid moves geodesically with rotation, 
then necessarily p = const. Eq. (9.4) implies then e = const, i.e. a geodesically moving 
and rotating perfect fluid whose rotation vector is covariantly constant must have constant 
matter density. Since A = may be assumed with no loss of generality (this leads only to 
redefining p), we shall assume this from now on. Then e = p and e + p = 2p is a conserved 
quantity. Hence, we may assume: 

n = e + p = 2p = const (9.5) 
in all formulae. In particular, (3.20) implies then: 

g = det(g a(3 ) = - hlA 2 = -(2p)- 2 . (9.6) 
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Using (9.6) and (3.23) in (9.1) we obtain, in the Plebariski coordinates: 

2 n (fa3,/3 - 9f33,a + 9 a p,z) = 0- (9.7) 

After a simple algebraic manipulation this set of equations yields the following result: 

533 = -A 2 = const, g a p jZ = 0, g 13jt = g 2 3,t = 0, 

923,x ~ 9l3,y = 0. (9.8) 

The second equation in (9.8) means that w a is a Killing vector, as should be expected from 
(9.1), (9.6) and (3.23). Eqs. (9.8) imply that 513 and 523 depend only on x and y, and that 
there exists a function T(x, y) such that: 

523 = T, y , 913 = F,x ■ (9.9) 

Since we assumed that rotation is nonzero, we know that 533 = —g a pw a w l3 /n 2 7^ 0, and so 
we are allowed to carry out the coordinate transformation: 

z = z'-F/g 33 , (9.10) 

that, in virtue of (9.9), will lead to: 

513 = 523 = (9.11) 
in the new coordinates. We have thus arrived at the metric form: 



ds 2 = (dt + ydx) 2 - h(t,x,y)[dx + k(t,x,y)dy] 2 - l(t,x,y)dy 2 - A 2 dz 2 , (9.12) 

where h, k and I are functions to be found from the Einstein equations and A is an arbitrary 
constant. 

From now on, the allowed coordinate transformations are (5.6) - (5.7), but with T = 
const. 

The components of the Einstein tensor will now be referred to the orthonormal tetrad 
implied by (9.12). The equation Gu = is integrated with the result: 

k, t = K{x, y)l^ 2 /h 3 / 2 - 1/h. (9.13) 
The equation G22 = k>P, with I elliminated by (9.6), is integrated with the result: 

h= [H 2 (x,y) + K 2 /(4k P )} 1/2 + H sm[2(K P ) 1/2 t + T{x,y)}, (9.14) 

where H(x,y), k(x,y) and r(x,y) are arbitrary functions. Now G\\ + G22 = 2kj> imposes 
an additional condition on (9.13) and (9.14) that leads to H = or: 

H = (A 2 - k/ P ) 1/2 K/{2k). (9.15) 

The case H = leads to the Godel solution (see below), so we shall consider the more 
general case (9.15). Then, from (9.14): 

h = [K/(2k)]{A + (A 2 - k/p) 1 ' 2 sm^npf'H + r]}. (9.16) 
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With such h, eq. (9.13) can be integrated with the result: 

k = [2A{K P ) l l 2 hY l {A 2 - n/p) 1 ' 2 cos{2{K P ) l l 2 t + t] + L{x,y), (9.17) 

where L is a new arbitrary function. The function I is then calculated from (9.6), and an 
explicit solution of Einstein's equations is determined by (9.16) and (9.17). 

The transformations (5.6) - (5.7) with T = const can now be used to simplify the metric 
tensor so that L = 0. The transformation that yields this is given in Appendix A. The still 
allowed coordinate transformations that preserve the property L = are given by (5.6) - 
(5.7) with T = const and with the additional condition: 

KF lX F, y +[k/(A 2 P K)]G, x G, y = 0. (9.18) 
With L = 0, from the equations Goi = G02 = one obtains further: 

r, x =2y( K p)^ 2 -Aip/n^K^, 

r, y =-( K /p) 1 / 2 K, x /(AK 2 ). (9.19) 
The integrability condition of (9.19) is: 

(p/K) 1 / 2 (AK), yy +(Kp) 1 / 2 [l/(AK)}, xx -2( K p) 1 / 2 = 0. (9.20) 

By the same method as was used in Ref. 4 it can now be shown that eq. (9.20) is at the 
same time the integrability condition for such a coordinate transformation (5.6) - (5.7) - 
(9.18) after which (see Appendix A again): 

K = ( K /A)y 2 , (9.21) 

and then (9.19) implies: 

T = c = const. (9.22) 

The value of c can be set arbitrarily by transformations of t of the form t = t' + const. To 
match Ref. 11 one should choose: 

c=-7r/2. (9.23) 
Finally, the functions h, k and I in (9.12) are thus: 

h = \y 2 {l + [1 - k/(M 2 )] 1/2 cob[2(k P )^% 

k = ^[(/ep) 1 / 2 /*]- 1 ^ - n/ipA 2 )] 1 ' 2 sin^M 1 / 2 *], 

I = (4p 2 A 2 h)~ 1 . (9.24) 

This is equivalent under a simple coordinate transformation to the Ozsvath class III solution 
from Ref. 11. 

The velocity field u a = 5q for this solution has nonzero shear. The shear will vanish if 
and only if: 

A 2 = K /p, (9.25) 
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and then the Godel solution in the form (5.11) results. 

The invariant definitions of the Ozsvath class III and of the Godel solutions given at 
the end of sec. 1 follow from the derivation in this section. 

10. Concluding remarks. 

These are the main results of the paper: 

1. With nonzero rotation, any Killing field, existing for a metric whose matter source 
inherits the symmetry, must have the form (5.5) when represented in the Plebahski coor- 
dinates. When <f>- a / 0, the Plebahski coordinates can be adapted to the Kiling field so 
that k a = 5 a 1 . 

2. When two of the generators of the group are spanned on the velocity and rotation 
vector fields, while the third one is not, the collection of solutions with a dust source is 
exhausted by two sets: 

a) The set of sec. 7, defined by a single differential equation (7.24), where the metric 
is (7.15) with A = k 13 = 0, F defined by (7.23) and y(Y) defined by (7.21). 

b) The set of sec. 8, where the metric is (7.15), with b = and the metric functions are 
defined by an autonomous set of first-order equations (8.1) - (8.3), (8.8) and (8.10) (the 
integral in (8.8) and (8.10) can be calculated if the variable is changed as in (8.9)). 

3. The solutions of Lanczos [1 and 19] and Godel [2] are limiting cases of both sets; 
their invariant definitions are given at the end of sec. 7 and of sec. 8. 

4. With no symmetries pre-assumed, if the source is a rotating geodesic perfect fluid 
whose rotation vector field is covariantly constant, then the solution of the Einstein equa- 
tions is the homogeneous (in four dimensions) Ozsvath class III solution [11]. If shear is 
zero in addition, then the Godel solution [2, see also 4] results. 

Note the modification that the results 2a and 2b introduce in theorem 3.1 of King and 
Ellis [42]. Those authors considered spatially homogeneous models in which the velocity 
field of matter was tilted (i.e. was not orthogonal) with respect to the hypersurfaces of 
homogeneity. Theorem 3.1 says, among other things, that there are no tilted models of 
type I and that tilted models of type II have zero vorticity. Evidently, this does not apply 
to the case where the hypersurfaces of homogeneity are timelike. The solutions of sec. 
7 are of Bianchi type II, they are "tilted" (because the velocity field is tangent to the 
hypersurfaces of homogeneity), yet rotation is not zero. The solutions of sec. 8 are tilted 
in the same sense, yet they are of Bianchi type I. 

Other solutions that have been published earlier will be mentioned where appropriate 
in papers 2 and 3. A general overview of literature on related subjects will be included in 
paper 3. 

The algebraic calculations for this paper were done with the help of the program Or- 
tocartan [43 - 44]. 

Appendix A. The transformation to L = in (9.17). 

A transformation of the class (5.6) - (5.7) with T = const changes the functions h and 
k to such ones that can be cast in the form (9.16) and (9.17), respectively, with the new 
functions K', t' and L' expressed through the old ones as follows: 

K' = K(F, X , +LG, X , f + kG, 2 x , /(A 2 P K), (Al) 
t' = t- 2(k P ) 1 ' 2 S + U, (A.2) 
L' = K'- l [KF, x , F, y , +2KLF, y , G, x > +KL + KL 2 G, X > G, y < +kG, x > G, y > /(A 2 pK)], (A3) 
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where S in (A. 2) is the function from (5.6) - (5.7) and U is determined by: 

cot U = 2A{kp) 1 ' 2 {2G, x , (F, x , +LG, x , )}' l [-K(F yxl +LG, X , ) 2 /{2k) + G, 2 X , /(2A 2 pK)}. 

(AA) 

Note that we are applying here (5.6) - (5.7) in reverse, i.e. with the roles of x a and x' a 
interchanged. The functions of the inverse transformation, denoted again by S, F, G and 
T, still obey (5.7). For consistency of all the formulae it is convenient to choose U from 
the segment (tt,2tt). Then, the limiting cases G, x > = and F, x > +LG, X ' = are included 
in (A. 4) as the limits U = 2ir and U = ir respectively (these limiting cases occur when 
L, xx = in the original coordinates). 

From (A. 3), the equation L' = turns out to be consistent with (5.6) - (5.7). In order 
to see this, one can solve (A. 3) and (5.7) for F, x > and F, y > and then impose the integrability 
condition F,^^ —F,^^ = 0. What comes out is a well-defined (though highly nonlinear) 
partial differential equation of second order for G whose coefficients depend only on K and 
L. 

The transformations preserving the property L = are (5.6) - (5.7) with (9.18), the 
latter easily follows from (A. 3). Equation (A.l) with L = then shows how K is changed 
by such a transformation; this is useful in showing that coordinates exist in which K, x = 
(see the remark after (9.20)). The proof is identical as in Appendix C to Ref. 4. Note that 
the conclusion in Ref. 4 is weaker than it could be: the K(x) (v(t) in Ref. 4) is determined 
up to an additive constant C. Hence, by a transformation of the form y = y'+ const and 
by an appropriate choice of C one can remove the linear and the constant terms in K (resp. 
v) so that K oc y 2 (resp. v oc t 2 in Ref. 4). 
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